Projection operator approach to lifetimes of electrons in metals 
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We present an alternative approach to the calculation of the lifetime of a single excited electron 
(hole) which interacts with the Fermi sea of electrons in a metal. The metal is modelled on the level 
of a Hamilton operator comprising a pertinent dispersion relation and scattering term. To determine 
the full relaxation dynamics we employ an adequate implementation of the time-convolutionless pro- 
jection operator method (TCL). This yields an analytic expression for the decay rate which allows 
for an intuitive interpretation in terms of scattering events. It may furthermore be efficiently evalu- 
ated by means of a Monte-Carlo integration scheme. As an example we investigate aluminium using, 
just for simplicity, a jellium-type model. This way we obtain data which are directly comparable to 
results from a self-energy formalism. Our approach applies to arbitrary temperatures. 

PACS numbers: 72.10.-d, 72.10.Bg, 72.15.Lh 



I. INTRODUCTION 

For several decades, the dynamics of excited electrons 
in metals has been the subject of intense research 
in theoretical and experimental solid state physics. 
Q, i, SIllSIESBISIiSl- These investigations are 
motivated by the fact that a lot of physical and chemical 
properties of metallic materials depend essentially on 
those dynamics [HI, [H, [iSl- Lifetimes of (photo)excited 
electrons in metals are always short (on the order of 
femtoseconds) but the immense progress in ultrafast 
laser technology now allows for an experimental deter- 
mination of such times, corresponding investigations are 
ongoing 

Today a number of methods are used to calculate life- 
times of electrons. Practically all of them are formulated 
within the framework of Green's functions (many-body 
theory) and aim at determining the self-energy, partic- 
ularly its imaginary part [l^ [3 [l3|- Many of them 
employ a screened interaction ("W") and a truncated 
expansion of the self energy in terms of this screened in- 
teraction ("GW-approximation") 0, i, [S E HO, HH, 
HI, m m m, m. III]. Xhe screened interaction is fre- 
quently obtained through a "random phase approxima- 
tion" (RPA) For a simple, sufficiently dense, ho- 
mogeneous gas of electrons interacting through coulomb 
repulsion (jellium model) an approach along the above 
scheme is even feasible analytically and yields a closed 
expression for the lifetimes close to the Fermi edge (see 
below) In a certain sense (which is described in 
more detail below p7|) ) this approach leads to lifetimes 
which may quantitatively be cornpared to experimental 
data on e.g., aluminium 0, [H, [23|. Of course timely 
state of the art approaches go beyond jellium and exploit 
not only the traditional self-energy formalism but also 
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density functional methods, etc. [23, |2J, |25|, |26l . 
Our approach is, in contrary, not based on Green's func- 
tions at all but on projection operator techniques. A 
main motivation of our work is to demonstrate that a per- 
tinent projective approach 28, 29, 30, 31, 32, 33, 34, 3j| 
is also capable to produce quantitative results on life- 
times. Our central formula from which the lifetimes are 
eventually calculated is in accord with expressions that 
may be derived within the above many-body approach 
(see below ([T5|) ). Furthermore it allows for an interpre- 
tation in terms of scattering events. This encourages a 
further development of projection techniques as alterna- 
tive quantitative tools for the investigation of relaxation 
and transport dynamics in condensed matter systems 
(find more on this at the end of Sec. lIVp . However, 
our approach starts form an effective model comprising 
pertinent quasi-particle dispersion relations and an ap- 
propriate screened interaction. The (generically subtle) 
provision of such a suitable effective model is not part of 
our present analysis, the effective model thus has to be 
supplied by other means. 

The article at hand is organized as follows: In Sec. HIl we 
give a very brief introduction to the time-convolutionless 
projection operator method [35] and apply it to a gen- 
eral interacting quantum gas thus obtaining an expres- 
sion for the electronic lifetime. In Sec. IIIII we evalu- 
ate this expression numerically for a "screened" jellium 
model tuned to describe aluminium. We compare our 
results to other available data and comment on comput- 
ing times. Eventually we close with discussion, summary 
and outlook. 



II. PROJECTIVE APPROACH TO 
OCCUPATION NUMBER DYNAMICS IN 
INTERACTING QUANTUM GASES 

To determine the lifetime of an electron initially occu- 
pying some momentum cigcnstate we analyze the dynam- 
ics of the corresponding occupation number. A formalism 
which allows for such an analysis is the TCL- method ^3^ . 
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In general the latter is a perturbative projection operator 
technique which produces autonomous equations of mo- 
tion for the variables of interest ( "relevant information" ) . 
The technique may be applied to quantum system with 
a Hamiltonian of the type H — Hq + XV where A has to 
be in some sense small SS^ In order to apply this method 
one first has to construct a suitable projection operator 
V. Formally, this is a linear map which projects any den- 
sity matrix g{t) to a matrix Pg{t) that is determined by 
a certain set of variables. These variables should match 
with the variables of interest. Moreover V has to fulfill 
the property of a projection operator, that is — V. 
For initial states with 'Pg{Q) = g{Q) the TCL scheme 
leads to a time-local differential equation for the dynam- 
ics of Vg: 

oo 

dtVg{t) = nt)Vg{t), T{t) = ^ X^'Tkit), (1) 

fc=i 

where the perturbative expansion used in the last equa- 
tions is in principle exact. However, for a description 
to leading order, which is typically and in our case the 
second order, one has to determine r2(i). Whether or 
not a leading order description will yield a reasonable re- 
sult is a somewhat subtle question [36|] but the expansion 
is well controlled and systematic, i.e, in principle higher 
order terms could be incorporated in a straightforward 
manner [ssl ] . A widely accepted indicator for the validity 
of the truncation is a clear timescale separation between 
the resulting relaxation dynamics and the decay of the 
correlation function, the latter being introduced below. 
However, here we are going to focus on the leading order 
and comment on the timescales below when we eventu- 
ally arrive at concrete lifetimes. In the literature [sH] one 
finds 

T2{t) = j dt'VC{t)L{t')V, (2) 



with C{t) — jr[V{t),^], where # denotes a placeholder 
for an operator which shall be inserted into the commu- 
tator. V(t) refers to a perturbation in the interaction 
picture. With ([1]) and ([2]) we obtain: 

t 

dtVgit) = j dt'VC{t)L{t')Vg{t). (3) 



Now for an concrete application we have to specify the 
underlying quantum model and a suitable projection op- 
erator. 

The systems we investigate are interacting quantum 
gases, here of the "spinless fermions"- type. The cor- 
responding Hamiltonians may be written as 

^ ekflk^k + ^(q)a[+qaLqai°^k, (4) 

k k,Lq 

Ho V 



where ek denotes a dispersion relation of free particles 
and V((\} the matrix elements of an interaction which 
depends on the concrete system. This Hamiltonian is 
of the above mentioned form as long as the interaction 
term V is in adequate sense "small" (see below). As will 
be demonstrated below (cf. Sec. IIII|) it is thus reason- 
able to choose for Sk pertinent quasi-particle dispersion 
relations and particularly for V^(q) an adequate screened 
interaction. Note that we neglect the spin quantum num- 
ber since the system we consider in this work is param- 
agnetic and without any magnetic fields the dispersion 
relation is the same for both spin channels. Below we 
are going to take care of this "spin degeneracy" in a 
very simple form (cf. text following (fT8|) . For the non- 
interacting many-particle system we may directly write 
down the wavenumber (momentum) dependent "single 
particle equilibrium density operator" as: 

qT ■■= T)ala; + (1 - /j(^, r))aja], (5) 

with /j(M,T) = (exp((e(j) - /x)/fcBT) + 1)"! being the 
Fermi distribution. Since we are interested in tempera- 
ture regimes close to T = OK but still T OK we can 
set the chemical potential ^ ~ ep- Further we abbre- 
viate fj{eF,T) as /j. The equilibrium density operator, 
again for the non-interacting case, of the total system, 
p*^', may be written as the tensor product of the single 
particle density operators, i.e. 

e'^'^ := (g) e- ^ note also g := gt' ■ (6) 

i i#j 

Here, for later reference, g denotes the total density op- 
erator of the system which does not contain the sub- 
space with respect to the momentum mode j, i.e. it is 
— We should, also for later reference, mention 

here that while g'^'^ is strictly speaking just the equilib- 
rium state of the non-interacting system, it is routinely 
considered to describe the single particle properties of the 
weakly interacting system more or less correctly. Thus, if 
single particle observables relax towards equilibrium due 
to the interactions (scattering), we expect them to relax 
towards values corresponding to g^"^. 
For the investigations of the dynamics of excited states 
we define an operator Aj (for the remainder of this paper 
|j) denotes the excited state) as 

Aj — (1 - - ./jajaj = ajaj - (7) 

which describes the deviation of the mode occupation 
number nj = ojaj from its thermal equilibrium. Now, in 
order to apply the TCL method to this model we con- 
struct a suitable projector as follows: 

Vgit)^g-'^ + ^TT{Aig{t)}g^A^, (8) 

with g{t) being the density operator which describes the 
actual state of the system, dj{t) :=Tr{Ajg(t)} denotes 
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the time dependent expectationvalue of Aj and a? := 
(1 — /j)^ + f? =Tr{A|}. It is straightforward to show 
that with the above definitions P is a projector and fulfiUs 
V^Q{t) = Vg{t). Note that Trj^.'^'Aj} = and Tr{Aje0 
Aj} = cr|. Before we eventually concretely apply ^ to 
our model we make the following approximation for an 
expression that appears in the computation of ([3]): 

C{t')VB{t) = '-[V{t'),e'''] + \[V{t')rQ® ^■a'^ 

« l[y(0,^®Aj]i|l, (9) 

The neglected commutator term essentially describes the 
dynamics of the equilibrium state of the non-interacting 
system. Eventually we are interested in a single particle 
observable. As already mentioned above, the equilibrium 
state of the non-interacting system is believed to reason- 
ably describe single particle observables in equilibrium 
even for weakly interacting systems. Since an equilib- 
rium state is constant, the above commutator should not 
significantly contribute to the relevant dynamics, thus we 
drop it. Keeping the term and performing all following 
steps yields eventually an expression which can explicitly 
shown to be indeed negligible in the weak coupling limit. 
For clarity and briefness we omit this calculation here. 
If we apply now ([2]) to ([H]) and make use of ^ we obtain 

t 

dtVoit) = \dtd;{t) • ^® Aj / dt'Q^'i 
t 

- / di'^Tr{[Aj[I/(i),[F(i'),^® Aj]]}^® Ajdj(i). 
^ 

(10) 

Multiplying both sides of pI7)) with Aj and taking the 
trace leads to: 

t 

MW = -^ / rft'Tr{[Aj[l/(t),[y(0,e®Aj]]}dj(t), 
^ 

rj(t) 

( 

where T-^{t) appears as a time-dependent damping rate 
of the mode j. If T^(t) turns out to be approximately 
time-independent, the usual exponential relaxation re- 
sults. With the substitution t' — t — t and exploiting 
(g) Aj , JJq] = [Aj , iJo] = as well as some trace proper- 
ties, we obtain for the rate: 

-t 

TjW = J^I rfr Tr{[y(T),g^Aj]-[T/(0),Aj]} , (12) 

^ C(r) 

where C (r) denotes the correlation function which is real 
due to the fact that both commutators are hermitian. 



The concrete evaluation of this expression with respect 
to our model is straightforward but somewhat lengthy. 
Thus the full computation is given in the Appendix, here 
we only give and discuss the results After exploiting the 
commutators within the trace we finally obtain for the 
rate: 

-t 

Ti{t,T) = l = -^J2 I rfr|y(q)pF(k,q,j,T) 
X cos((wk-i-q + Wj_q-a;k-Wj)r) 

= ^ E q. j: T)-t- sine (w(k, q, j)0 , 

J k,q 

(13) 

with: 

w(k,q,j) := t^k+q + t^j-^q-^k- Wj, 
F(k,q,j,r) := (l-/j)(l-/j-q)(l-/k+q)/k 

i=^l(k,q,j,T) 

+ ./j/j-q,/k+q(l - /k) 
F2(k,qJ,T) 

(14) 

and sinc(a;i) denotes the sinus cardinalis. Obviously the 
integral t-smc{LLit)dijj is independent of t. Further- 
more the function gets more and more peaked with in- 
creasing t such that, as wellknown, 

sinM^^^(^), (15) 

t — ^oo (xJ 

Hence, since dispersion relations are smooth functions of 
the wavenumber, we expect the rate T-j(t,T) to become 
indeed time-independent for times larger than Tc, if I/tc 
is an energyscale on which dispersion relations may be 
linearized. Thus for times t larger than Tc we may with 
good precision approximate (here we neglect the factor 
since for temperatures T « OK it is w 1): 



^S(T) = 1 = ^ ^ |y(q)p5 + ej_q - ek - ej) 

j k,q 

X {(1 - /j)(l - /j-q)(l - /k+q)/k + /j/j-q/k+q(l - /k)}- 

(16) 

This expression is one of our main results. In princi- 
ple it allows for a direct calculation of lifetimes for any 
fermionic system with given quasi-particle dispersion re- 
lations and screened scattering term. Very similar formu- 
las can be found in textbooks in the context of transport 
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and relaxation, see, e.g., [lO|, [llj. They are often de- 
rived on the basis of an ad hoc apphcation of Fermi's 
golden rule. A closer look reveals that such an expres- 
sion can also be obtained from an analysis along the lines 
described in the introduction ("RPA-GW") by using the 
static-limit form of the screened interaction. Since our 
further quantitative determination only consists in a nu- 
merical evaluation of (|16p the outcome is equivalent to 
the one obtained by the above treatment. And since, as 
outlined in the following, (|16p is in accord with a stan- 
dard scattering interpretation, obviously both, the pro- 
jective and the above version of the many-body approach 
amount more or less to the counting of scattering events. 
The contributions to the decay rate corresponding to Fi 
and F2 allow for an intuitive interpretation, at least for 
low temperatures. 

Fi : This term accounts for the decay of an electron from 
a momentum mode j above the Fermi sea, cf. Fig. ([T^). 
Due to the factor (1 — /j) it only significantly contributes 
to the occupation number dynamics of such modes that 
are unoccupied in equilibrium. Those occupation num- 
bers may only deviate from equilibrium towards an ex- 
cess of electrons. According to the other three factors 
only those summands contribute that correspond to the 
electron at j colliding with an electron from within the 
Fermi sea k, such that the post-collision momenta k + q, 
j — q lay in the unoccupied region above the Fermi sea. 
F2: This term accounts for the decay of a hole from a 
momentum mode j within the Fermi sea, cf. Fig. IHd). 
Due to the factor /j it only significantly contributes to 
the occupation number dynamics of such modes that are 
occupied in equilibrium. Those occupation numbers may 
only deviate from equilibrium towards a shortage of elec- 
trons, i.e., holes. According to the other three factors 
only those summands contribute that correspond to two 
electrons from within the Fermi sea k + q, j — q colliding 
such that one post-collision momentum k lays in the un- 
occupied region above the Fermi sea and the other lays 
exactly at j such as to fill up the hole. 



a) 




FIG. 1: Schematic representation of the underlying collision 
processes in momentum space as described in the text (kp 
denotes the Fermi momentum). 

a) A collision process through which an excited electron at 
momentum j vanishes from its initial momentum mode. 

b) A collision process "filling" a hole within the Fermi sphere 
at j. 

The dashed circles denote the possible outgoing momenta un- 
der momentum and energy conservation. 



nian of the model is given by: 



III. APPLICATION TO A JELLIUM MODEL 
WITH SCREENED INTERACTION 



In this section we now apply our result for the de- 
cay rates to a jellium modell featuring a Thomas-Fermi 
screened interaction. The latter will eventually be tuned 
to correspond to aluminium. However, to repeat, the 
main intention of this work is not to calculate decay rates 
in aluminium with extreme precision, but to concretely 
demonstrate the feasibility of our method. The Hamilto- 



IE 



aL„a|_„aiak, 



(17) 



where Vt denotes the volume of the solid and q^p the 
so called Thomas-Fermi wavenumber which is related to 
the Fermi wavevector and the Wigner-Seitz radius by 
{qTF/kpf' = 0.665r5 (with = (4;^^)^^, oq being the 
Bohr radius, ep = (9^/4)i-^[ryd], kp = (9^/4) 
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qTF = (I^/tt) 3 ^--i=). Note that our model only com- 
prises electrons, no phonons. Thus the result on the de- 
cay rate has to be compared to that part of the total 
decay rate that stems from electron-electron scattering 
only. In real aluminium there is evidence that the total 
lifetime is also significantly shortened due to electron- 
phonon scattering 0, We apply now (fTC|) to this 

model which yields: 



1 47r v-^ 



An ■ 
li 



X 5 (ek+q 



k,q 



J^(k,q,j,T) 

(18) 



where the auxiliary factor 2 arises from the fact that for 
each k we have two one-electron states (one for each spin) 
which is taken into account by this additional factor. If 
we now replace the sums in ([TS]) by integrals by the rule 
^/(k) ^ (27r)~3f7/dk/(k) we obtain for T^{T): 



rj(r) 



5 (^k+q 



dkdq 



dkdq 



i^k - e-i) 

i^(k,qj,r) 



i^(k,qJ,T) 



WhelTT^JJ ^(q2+q2,^) 



(5 (£k-|-q 



'j-q 



- £k - £j) ■ 



(19) 



For the numerical calculation it is advantageous to trans- 
form the momenta to dimensionless parameter thus we 
introduce coordinates relative to the Fermi momentum: 
h ^ kp ■ k', thus q — > /cf • q' and thus j ^ kp ■ j'. 
Applying all these substitutions to (fT9|) leads to: 

=4 



rj'(r) = 



dk'dq' 



j"(k^q^j^^) 

(q'2 -I- 0.665rs)2 



X <5((k' + q')^ + (j'-q')^-k'2-r), (20) 

with /k' = (exp (/3eF(k'2 - 1)) + 1)"^ and = 9.1 • 
10^'^^ kg the free electron mass. 

For the numerical evaluation of the last expression we 
approximate the delta distribution by a suitable non- 
singular, e.g., Gaussian-type function: 

<5,(w)« ^=e-i^, (21) 



crV27r 

where a denotes the standard deviation which has to be 
in some sense small. More details on this somewhat sub- 
tle approximation are given below. Thus we get the fol- 
lowing integral for the rate: 



dk'dq 



, F(k',q',j',r) 



(q'2 -t- 0.665rs)2 



1 



exp -^((k' + q')2 + (j'-q')^-k'2-j'^ 




FIG. 2: Comparison of the lifetimes of excited electrons in 
aluminium as arising from electron-electron scattering only. 
Displayed is the regime close to the Fermi edge. The open tri- 
angles correspond to data obtained from from the approach at 
hand, i.e., integration of (|22p (T= lOK). The solid circles are 
experimental data corrected for transport effects taken from 
([j) while the solid line is the analytical result from Solid 
diamonds are the results of an DFT-GW calculation as shown 
in The number of sample points for the Monte-Carlo in- 
tegration of (|22|) is N=2 ■ 10^, the broadening parameter is 
o-=l/25. 



For aluminium we choose rg — 2.07. The six dimen- 
sional integrals are solved numerically without any fur- 
ther simplification using a standard Monte-Carlo package 
as implemented in the Mathematica code. Of course this 
specific integral could be evaluated in other ways, how- 
ever, to demonstrate the feasibility of our approach in 
general we proceed as indicated. 

As one can see in Fig. [5] there is rather good agreement 
between our results, other theoretical approaches and ex- 
periment. Our data is denoted by open triangles. The 
solid line corresponds to the many-body approach based 
on jellium [l| as outlined in the introduction. The solid 
diamonds denote the result of a more sophisticated many- 
body approach which takes the lattice into account and 
exploits density functional theory Q . The solid circles in- 
dicate the parts of the measured decay rates that are at- 
tributed to direct electron-electron scattering, i. e., after 
removal of transport effects according to [1]. FigO shows 
the analytic result from Q (263rs^^^(£-eF)"^[eV]2[fs]), 
which is supposed to be valid close to the Fermi edge, 
boldly continued to all energies (solid line). Further- 
more results of our approach for all energies are displayed 
(dots). Obviously there are deviations for electrons at 
higher energies while the agreement remains very good 
in the limit of "low-energy holes" . 

However, a comment should be added here. For this 
more or less realistic model we get lifetimes on the or- 
der of some femtoseconds. The decay time of the cor- 
relation function p2p is, very roughly, on the order of 



(22) h/e^ 



with 



being the bandwith. For about lOeV 
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T(e') 




0.0 0.5 1.0 1.5 2.0 



FIG. 3: Comparison of the logarithmic lifetimes of excited 
electrons (above e' = 1) and holes (below e' = 1) in alu- 
minium as arising from electron-electron scattering only. Dis- 
played is a wide regime around the Fermi edge. Data are 
plotted over rescaled energy, e' — e/ep. Displayed are re- 
sults obtained from Fermi-liquid theory as cited in ([J]) (solid 
line, T=OK) and from numerical integration of (|22|l (dots, 
T=10K). The number of sample points for the Monte-Carlo 
integration of (HU is N= lO'' and a = 1/10. 

this yields ca. half a femtosecond. Thus the separation of 
those timescales, which has been mentioned in Sec. Illll as 
a criterion for the truncation performed above, is not as 
clear as often in other fields, such as, e.g., quantum op- 
tics. This indicates that such models, at short lifetimes, 
are barely in the Markovian, weak coupling regime and 
hence memory effects and/or higher orders may have sig- 
nificant influence. 

To the choice of cr: Obviously a smaller a leads to a 
better approximation of the ^-function which should be 
the correct weight distribution at least in the long time 
limit. However, recall the above discussion of the time- 
independence of the decay rate. For analogous reasons 
larger a should leave the result unaltered, as long as a 
remains small enough to allow for a linearization of the 
dispersion relations on the scale of cr. A large a is nu- 
merically favorable since the larger a is, the larger will be 
the fraction of the Monte Carlo points that significantly 
contribute to the integral. And of course this yields a 
decreasing statistical error. Thus, for a given statistical 
integration error, a larger a simply implys a longer com- 
puting time. Hence finding the best a is an optimization 
process that should be done carefully. However, to name 
a number, the computation time for one of the lifetimes 
as displayed in Figs. HE] is about an hour. 

IV. SUMMARY, CONCLUSION AND 
OUTLOOK 

In this paper we considered the lifetimes of (quasi- 
)particles or holes in interacting quantum gases (only 
electronic part), using a projection operator technique. 
This yields a formula for the decay rates into which es- 
sentially the pertinent, effective quasi-particle dispersion 



relations of the particles and their screened interactions 
enter. This formula turns out to be in accord with an 
expression that may be found from a certain implemen- 
tation of the self-energy formalism . The rates are even- 
tually given in terms of integrals which can be cast into 
a form which is well suited for a Monte Carlo integra- 
tion scheme. While this work essentially aims at demon- 
strating the feasibility of this approach in general, the 
method has been concretely applied to a jellium model 
featuring a Thomas- Fermi screened interaction (tuned for 
aluminium) as a simple example. Here it yields reason- 
able results while requiring moderate computational ef- 
fort. This motivates an application of the approach to 
more complex systems. However, the results on life and 
correlation times indicate that such systems are, for short 
lifetimes (high electron energies, etc), barely Markovian 
and thus the decay may not even be strictly exponential. 
This hints at a necessity to include higher order terms in 
future investigations in this regime. 
The approach at hand aimed at generating an au- 
tonomous, linear equation of motion for a single electron 
occupation number (jlip . However a slight modification 
of the projection used here may directly yield linear equa- 
tion of motion for all electron occupation numbers, i.e., a 
linearized Boltzmann equation. As wellknown, the latter 
is a traditional starting point to investigate, e.g., trans- 
port properties. To those ends one would use a projection 
very much like the one discussed here ([8]) but summed 
over all occupation numbers j. The reasonable results 
on lifetimes presented in this work may be viewed to en- 
courage further investigations in that direction. 



APPENDIX A 



In this section we show the derivation of p6p . The 
main work is to exploit the two commutators and finally 
the trace. First we exploit the commutator [F(0), Aj]: 



[f(0),Aj] = 2 E^(q)[4+q«iVi«k,al«J-/j] 

k,l,q 

= ^ 51 ^(^)[°k+q"^-q°^l"k,a]aj]• 
ka,q 



Since the commutator is zero for j 7^ k + q, 1 — q, k, 1 
we just have to regard cases where one of the indices is 
equal to j and note that OittjOi = Oi, ajCjOi = 0. From 
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this follows that: 

k,q 

k,q 

+ ^ 13 ^(q)aj ajai-qaj+q. (Al) 

^ l,q 

With suitable index shifts and the fermionic commutator 
relations we finally obtain for the commutator 

[F(0), Aj] = ^ T/(q) (al+qa]_qakaj - ajaj,aj_qak+q) • 



where the annihilation and creation operators act 
only on the respective subspaces of the tensorprod- 
uct of the single density operators. With a\a\o\ = 
0, atfliaj = aj, the rules above and u(t) := 
exp (£k+q + ei-q - £k - £i) t) it follows: 



k,q 



(A2) 



Now we deal with the second commutator [V(t), q® Aj] 
where we first regard the case j ^ k + q, 1 — q, k, 1 (we 
abbreviate := 1 — /i): 

[l/(T),e® Aj]^J 

= \Y^ ^(q)k+qW«Lq(^)ai(^)ak(r), ~Qs ® Aj] 

k,l,q 



= i ^ y(q)eH(^fc+<i+=i-<i-=k-£i)T 

^ k,l,q 

>^ ["k+q^l-q'^l^k, (/iajfli + Siaifli) ® Aj], 



[F(T),e® Aj]^J 

^ X] ^(q)"('^) (/k/l5l-q5k+q - 5kffl/l-q/k+q) 



k,l,q 

X ak+qaLqOlOk 



Mj,k, 
l,k+q, 
1-q 



(A3) 



For the cases where one of the indices k + q, 1 — q, k, 1 
is equal j we obtain analogous: 



[F(r),^® Aj]=J = ^^V(q)a^q(r)a]_^(T)ak(T)aj(T) (g) ^f' (5j5j-q5k+q/k + /j/j-q/k+q5k) 

k,q if^j.k, 

j-q> 

k+q 

+ \^y{^a\+^^'^)a\_^{T)a-AT)ax{T) ^f" (a5l-q9j+q/l + /j/j+q/l-qfll) 

l,q iT^j.l, 

1-q, 
j+q 

- '^^(q)«k+q('^)4('^)'^k(T)aj+q(T) (g) (ft+q^k^j /k+q + /j/k/j+q5k+q) 

k,q i#j,j+q, 

k,k+q 

- \^'^{'^)a\{T)a\_^{T)a-;^-q(T)ax(T) (5j-qffl5j/l-q + /j/l/j-qffl-q) • 



l,q 



iT^jJ-q. 

1,1-q 



Again with suitable indexshifts and substitutions it fol- 
lows for the commutator: 



(A4) 
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[V{t), ~e ® Aj]=J = ^ Via) {al^^{T)a\_^{T)a^.{T)as{T) - a](r)aj,(r)aj_q(T)ak+q(r)) 



k,q 



X (aft-qSk+q/k + Zj/j-q/k+qSk) (X) 



j-q, 

k+q 



(A5) 



Now we exploit the trace: 

C(r) = Tr {[V(r), ^ ® Aj][l/(0), Aj]} - Tr { {[V{t), g ® Aj]^J + [V{t), g ® Aj]=j) [y(0), Aj]} 
= Tr{[y(r),^®Aj]^J[f(0),Aj]}+Tr{[l/(r),e®Aj]=J[F(0),Aj]}, 



A(r) S(r) 

where we split it into two parts and exploit them respectively. 



k,l,q 



4+q«i-q«i«k Pi'^ «) Aj (4_,_yoj_yaxaj - a]4aj_yax+y 



Mj,k, 
l,k+q, 
1-q 



= JX^nq)nyMT)G(k,l,q,T)Tr< 



k,l,q 



,69 K?, A.^t „t 



- J E ^(q)^(y)«(T)G(k, 1, q, T) Tr < 



k,l,q 




(A6) 



with G'(k, 1, q, T) = /k/i.9i-q.9k+q — 5'kfl'i/i-q/k+q- We focus now on the traces. For taking the trace we use the 
occupation number representation. 



Tr ( 



^k+q"l- 



l,k+q, 
1-q 



k+q"'l-q"l°'k 



■ ^ (g) Aj • a],+ya]_yaxaj|ni, ...,nr} 



(A7) 



Since under consideration no one of the indices k, 1, k + q, 1 — q is equal j this trace is zero for cases with y 7^ which 
is valid for II also. The case y = must be analyzed independent. 
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I for y = (now we write down the sum again) we have: 

\ E V{ci)V{0)u{t)G{\^, 1, q, r)Tr [al^^a\_^a,a^ • S ® Aj • 4a]a,,aj} . 



(A8) 



k,l 
q,x 



Here we have two summands: the case where k + q = 1 and l + q = l^q = 0: 

\ E F(q)T/(0)«(T)G(k, 1, q, T)Tr [al^^a\_^a,a^ • S ® Aj • ala^a^a^] 



k,l 
q,x 



\ ^(q) V^(0)u(T)G(k, k + q, q, T)Tr {ai,+qal,ak+qak • S ® Aj • aia\a^ai } 

k,q,x 

i ^ l/(0)F(0)w(T)G(k,l,0,T)Tr{a[aJaiak-S® Aj -atajaxaj} =0, 



(A9) 

since G(k, k + q, q, T) = G(k, 1, 0, T) = 0. For II the argumentation is analogous. Thus there is left just one more 
possibility: the case j = x + y for which we obtain from (jA6[) : 



k,l,q,x 

SO that finally follows that A{t) = 0. 



\ E ^(q)^(y)«(T)G(k,l,q,T)Tr{a^qaLqa,ak-S0Aj.a]a]_yaj_yaj} 

k,l,q,x 

\ E nq)^(y)^MG(k,l,q,T)Tr{a^qaJ_qa,ak-5® Aj -ajaj^yaj^yaj} =0, 



(AlO) 



For -B(t) we have: 



B{t) = ^F(q)F(y)^^(k,qJ,T)TrNa^q(r)a]_q(T)ak(T)aj(T)-a](r)ai,(T)aj_q(r)ak+q(r)) (g) 



k,q 



j-q, 

k+q 



-2 ^ |l^(q)P (.gjffj-q.gk+q/k + /j/j-q/k+q.9k) COS ((wfe+q + CJj-q - Wfc - t^j)T) , 
k,q 

I 



(All) 



from this follows (fT6|) . 
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